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Several recent experimental studies have confirmed the possibility of switching the magnetization 
direction in the small magnetic domains by pumping large spin-polarized currents through them. 
On the basis of equations proposed by J.Slonczewski for the case when magnetization of the do- 
mains is almost uniform, we analyze the stability and switching for two types of magnetic and shape 
anisotropies of a magnetic domain in a nanowire and find qualitatively different behavior, including 
different shapes of bistable regions. Our study is analytic as opposed to recent numeric work. As- 
sumed anisotropies can be realized in experiments and our predictions can be used to experimentally 
test the theory of spin-transfer torques. Such test would be especially interesting since alternative 
approaches are discussed in the literature. 



If the size of the pieces is larger then the domain wall 
width, the magnetization may be not uniform throughout 
the piece. In this case equations from || must be used in- 
side each piece to determine the magnetic configuration. 
In this work we will assume that magnetic pieces are suf- 
ficiently small and treat them as magnetically uniform. 
Then the total anisotropy, including magnetocrystalline 
and shape contributions, will be given by a symmetric 
tensor K\£\ 

The problem simplifies if one of the magnetic pieces is 
much longer then the other. At a given current density j 
the spin-transfer torque T st is proportional to the cross- 
section of the piece while the torques T a due to anisotropy 
are proportional to its volume. The ratio T st /T a ~ 1/L, 
where L is the length of the piece, thus the small piece 
will be affected by the spin-transfer torque starting from 
a much smaller value of j. We can therefore neglect the 
effect of T st on the large piece, called a polarizer, and 
assume its magnetization to be constant. Magnetizations 
of the large and small pieces will be denoted as M^ and 
Mg respectively. 

To test the theory, one would like to be able to control 
the direction of the magnetization of the polarizer with 
respect to the anisotropy directions of the small piece. 
The easiest way to change M^ is an application of an 
external magnetic field B. Of course, B will also act on 
the small piece and must be taken into account in the 
equations of motion. The properties of a system with a 
given anisotropy K^ can ultimately be presented as a 
phase diagram in the (j, B) space with spin-transfer ef- 
fects determined by the magnitude of the current and by 
Mi = Mi(B). In this paper we calculate a section of 
such diagram for certain directions of B and anisotropies. 
Note, that for technical applications in the memory writ- 
ing process one is interested in finding anisotropy tensors 
which satisfy two conditions: (a) there is a section of the 
phase diagram at a fixed external field B = Bo where Ms 
is bistable at j = 0; (b) by passing a current j ^ one 



Recently considerable experimental interest [|l|-£§ has 
been shown in the torques created by spin-polarized cur- 
rents in a magnet. This interest is fueled in part by the 
proposals of developing a convenient writing process for 
advanced metallic magnetic RAM || where the reading 
process will be based on the magnetoresitance or other 
effect || . A general theoretical framework for the de- 
scription of such "spin-transfer" torques is set in |Q-Q] . 

One of the particular experimental setups in which this 
effect can be studied is a thin (< 300 nm) normal metal 
wire with two magnetic pieces embedded in it (see Fig.l). 
If the distance between the magnetic pieces does not ex- 
ceed the spin diffusion length in the normal spacer be- 
tween them, a current passing through the wire will in- 
duce spin-transfer torques in both magnets. Such setup 
was originally considered in M for the case when both 
magnetic pieces are isotropic and their magnetizations 
are initially not collincar. It was predicted that both 
magnetizations will perform rotation in a fixed plane 
keeping the angle between them constant. 

B ^ M, 
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FIG. 1. Model experimental setup. Current j is passed 
through the nanowire with two magnetic pieces (shaded areas) 
with magnetizations M_l and Ms. External magnetic field B 
can be applied in an arbitrary direction. 

However for a real material one must also take into 
account magnetocrystalline and shape anisotropy. In 
this paper we present the first exact solutions with 
anisotropies taken into account. Attempts to incorpo- 
rate them were made in the experimental reports pfj] 
and an extensive treatment of a particular experimental 
situation was given in a numerical study [hoi . 



can switch between these two metastable states. However 
for the purposes of testing the theory it is better to start 
with the cases where the phase diagram can be calcu- 
lated exactly and compare theoretical and experimental 
results. 

In p0| | stability of the zero-current equilibrium position 
of Ms was studied by an approximate analytical method 
and by numerical solution of a differential equation for 
the time evolution of Mj. The latter gave the final po- 
sition of magnetization and the details of the switching 
process. Due to the approximation used, there was a dis- 
crepancy between the switching thresholds predicted by 
two methods. In our study we find the equilibrium posi- 
tions of Mg for each (j, B) point and analyze their sta- 
bility exactly. Knowing the nature of all equilibria we are 
able to construct the topology of the time-evolution flow 
of Ms and predict qualitatively the overall behavior of 
the system, including existence of the stable cycles. Such 
stable cycles were predicted in |13| and observed numer- 
ically in |Hj. Due to energy dissipation they would be 
impossible without the current. However for j ' ^ there 
is a constant supply of energy which feeds the periodic 
motion of Ms- 

Dynamic Equation for the Small Piece The mag- 
netic energy of the small piece is a sum of the intrinsic 
anisotropy term, shape anisotropy term, interaction with 
external magnetic field and exchange interaction with the 
large piece. Approximating the shape of a monodomain 
small piece by an ellipsoid we can write 



y = ~ (-K[l ntr) n t n k + AnM t N lk M k ) 
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where Nik is the demagnetization tensor pf , Mi = M$i - 
magnetization of the small piece, K^ mtr ^ - magnetocrys- 
tallinc anisotropy tensor, and J ex is the exchange cou- 
pling between the pieces. Vectors n and s are unit vec- 
tors along the magnetization of the small and large pieces 
respectively. According to M, the modified Landau- 
Lifshitz equation for M has the form: 



+ 2^ 5 (P, S .n)[n,[ S ,nI 
V 2 e 



(1) 



a n, n 



where 7 = g^s/h is the gyromagnetic ratio, V and A 
are the volume and cross-section area of the piece, the 
last term represents Gilbert damping, P is the degree of 
spin-polarization of the electrons coming out of the large 
piece and g(P, s • n) is the function derived in M 
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Equation (nj) can be rewritten in terms of n 



with rescaled coefficients 
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(4) 



The behavior of the small piece will be completely deter- 
mined by the directions of s and h with respect to the 
principal axis of the anisotropy tensor K . Dependence 
s = s(h) is given by the properties of the polarizer. We 
will assume that the damping is small and expand the 
solutions in a. 

It was found recently |l|] that cobalt nanowires grow 
with intrinsic easy axis perpendicular to the wire for large 
wire diameters d > 50nm and with easy axis along the 
wire for smaller d. As for the shape anisotropy contribu- 
tion, it will be an easy axis along the wire if the length 
of the small piece Ls » d and an easy plane perpen- 
dicular to the wire in the opposite case Ls << d. In 
P,flOf a wire with d ~ 100 nm was used. This compli- 
cates the situation because for any Ls all three principal 
components of the total anisotropy tensor are different. 
On the contrary, for d < 50nm wires anisotropy is al- 
ways uniaxial. For Ls » d it is an easy axis along 
the wire. For Ls « d the total constant is given by 
K = /-sT( iTltr ) - 4:nM 2 . If M is sufficiently large to ensure 
K < 0, one has an easy plane anisotropy. This is the case 
for cobalt K = 5 • 10 6 erg/sm 3 and M = 1.4 • 10 3 emu. 

We are going to consider two experimental situations 
with thin wires. 

Axial Case. Assume that the polarizer is character- 
ized by an easy axis anisotropy along the wire. The small 
piece has uniaxial anisotropy with respect to the same 
axis, but the anisotropy constant K can have either sign. 
Next, assume that the external magnetic field is also di- 
rected along the wire which itself is oriented along the 
z-direction. Such situation with K > was considered 
before in |l|] using a different method. 

First we rewrite vector equation (|3|) in terms of the 
polar angles of n: </> and 9. That gives a system: 
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To find equilibrium positions one must solve 



v^ = sin 6(h + K cos 9) = 
vq = Ig(P, cos 9) sin 9 = 



(5) 



(6) 



[h, n]i + K ik n k +Ig(P,s ■ n)][n, [s,n]]j +a[n,h] (3) 



When / ^ the only stable positions of n on the unit 
sphere are the North and South poles, independent of the 
magnitude of current. 

To determine the stability of equilibria we linearize the 
r.h.s. of (||) in small deviations. At 9 = (0, 7r) one has to 
use local non-singular coordinates x = 9 cos <f>,y = 9 sin (f>. 
Linearized equations have the form 
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The nature of equilibria is determined by the eigenvalues 
fj, of the dynamic matrix D. Real /i's imply a center or a 
saddle, while for complex-conjugate eigenvalues the sta- 
tionary point is a stable focus for Re/i < or an unstable 
one otherwise. 

Without showing the algebra we write down the result. 
For the North pole 

fi N = -Ig(P, l)-a(h + K)±i\h+ K - aIg{P, 1)| (8) 

and therefore it is a stable equilibrium for 

a(h + K) 



For the South pole we get a stability criteria 
a(h-K) 



j < 



g(P, -i) 



(9) 



(10) 



The regions of stability of the North and Sough poles are 
shown on Fig. 2. It is important that there is a region on 
the diagram where both equilibrium points are unstable. 
This necessarily means that there exists a stable cycle, 
around which n performs a periodic motion. 




FIG. 2. Switching phase diagram for the axial case. (A) 
easy axis in the small piece; (B) easy plane in the small piece. 
The left boundary (vertical dashed line) represents the switch- 
ing field of the polarizer. Regions are marked as: N - only 
North pole is stable, S - only South pole is stable, N/S - 
bistable region, P - precession region. Time-evolution flow 
on the projected sphere (in which N is mapped to infinity 
and S to the origin) is shown on the insets. 

In general, stability analysis can not give any infor- 
mation about the shape of the cycle, but in the ax- 
ial case it is easy to anticipate that it will be given 
by 9 = 9 a . From (||) the value of 9 is given by 
I g(P, cos 9 ) + a{h + Kcos9o) = and the precession 
frequency by ui — \Ig(P, cos 9 )/a\. When h is de- 
creased through the precession region at fixed current, 
#o changes continuously between and 7r and u> increases 
from uj n = \Ig(P,l)/a\ to cj s = \Ig(P,—l)/a\. For the 



K > case / > 2aK/(g(P,-l) - g(P,l)) in th e pre- 
cession region and the frequency has a minimal value 
2Kg(P, l)/{g(P, -1) - g(P, 1)) of the order of FM reso- 
nance gap. For K < the frequency can go down to zero 
since there are 1 = points in the precession region. 

Planar Case. We now consider a situation when the 
polarizer is characterized by an easy-plane anisotropy. 
The small piece has no choice but to have K < too, 
since by definition L$ « Ll and shape anisotropy must 
be even more important for it. In this experiment we 
direct the magnetic field perpendicular to the wire. As a 
result, we always have s along h and there is no meaning 
to negative values of h. We choose s and h to point in 
the x-direction. 

Functions v^ and vg for the planar case differ from the 
axial one. Working out their form we get a system of 
equations for the equilibrium points. 



v<j) = K cos 9 sin 9 — h cos 9 cos <f> — Ig • sin <j> — 

vg = h sin <fi — Ig ■ cos 9 cos <\> = 



where g = g(P,sm9cos(j)). It is that dependence of g on 
both 9 and </> that brings the main complication into the 
solution for the planar case. 

Both equations are satisfied for cos 9 = sin0 = 0. Con- 
sequently there is a pair of equilibria: point A (</> = 0, 9 = 
ir/2) and point B (<f> = n,9 = n/2), positions of which do 
not depend on the current magnitude. There are however 
additional equilibria given by 



K sin 9- cos <f>{h 2 + I 2 g 2 )/h 
h sin <j) — Ig- cos </> cos 9 
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To deal with the 
change 



, 9 dependence of g we make a variable 



cos <p sm o = x 
sin = y 



(12) 



Then g = g(P,x). It is an important property of the 
model and the planar case geometry, that this substi- 
tution succeeds in reducing system ( pd| ) to a relatively 
simple decoupled equation on x: 



1 + (vxg(P, x)f 



x 



(13) 



with xq = h/K < and r/ = I /h. 

Investigation of (O) shows that as r\ increases, there 
will be one, two or zero roots. We will see below, that 
switching happens already at small currents, where only 
one root exists and where we can expand the solution 
as x — xq + r\x\ + . . .. This root defines two symmetric 
equilibrium points C\.i located as shown on Fig. 3 inset. 
Their positions depend on the magnetic field. As h grows, 
both C points move gradually towards B and merge with 
it at h = \K \. At higher r\ we observe additional roots - 
the full (/ — h) diagram will be published elsewhere. 



Stability analysis shows, that point A becomes unsta- 
ble at a negative current with a magnitude I sw ~ a, 
similar to the axial case. That justifies the expansion of 
x above. As in the axial case, we only present the eigen- 
values of the dynamic matrix at all equilibrium points. 

For point A 

HA = -Ig(P, 1) - a (~ K 2 + 2k) ± i^/h(-K + h) (14) 

Therefore A is a focus point. Its stability is given by 

-K + 2h 



Linearization near point B gives 
2ft + if 



(15) 



H B =Ig{P,-l) + a- 



± ^-h(h + K) (16) 



Therefore point B is a saddle point for h < \K\ and is a 
focus for greater field. In the latter case B is stable if 



I < —a 



2ft + K 



MP, -1) 
Treatment of both points C% gives identical results 



(17) 
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1 + COS 2 I 
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2aK + oi- 
ls. 



± iy/K? - ft 2 + 0(1) 



In the region of its existence (ft < \K\) point C is 
always a focus. The stability region is given by 



I < -aK 



X (2-x^) 



2 + fg(P,x )(l-xo) 2 x 



(18) 




Inequalities ( |15| , pL7|Jlq ) determine the switching phase 
diagram in the small a limit (Fig. 3). In the planar case 
only bistable regions exist and there is no precession re- 
gion on the phase diagram. We also note that the basins 
of attraction of C\ and Ci are increasingly interpenetrat- 
ing near the boundary separating them from the basin of 
attraction of A (Fig. 3). In switching from A to G the 
choice between C\ and C2 is random. 

Discussion. Switching patterns depend crucially on 
the magnetic anisotropy and the direction of polariza- 
tion of incoming current. Our predictions for the axial 
and planar cases can be used to experimentally test the 
theory , in particular the accuracy of the factor g(P,s ■ n) 
M. This is especially interesting because alternative de- 
scriptions of current driven excitations are put forth in 
the literature [ |14|Jl5| . Since currently magnetization di- 
rections are not experimentally measured but rather in- 
ferred from the resistive state of the wire, it is important 
that there are qualitative differences between axial and 
planar cases. In the axial case all boundaries on the 
(/ — h) diagram are straight lines, while in the planar 
case one of them is curved. Current sweep through the 
precession region in the axial case will show three resis- 
tive states without hysteresis, while a sweep through the 
bistable region gives two jumps with a hysteresis. If Mj 
is measured directly, one will see that it rotates by 180° 
degrees in the axial case and by a magnetic field depen- 
dent angle ir/2 + 80(h) in the planar case. The precession 
state is a good candidate for observation with a magnetic 
force microscope. 

The switching current density can be calculated by 
substituting (g) into (g [[(]) and (|H||^,|l8|) . As a charac- 
teristic value one can take 



\K\ 



g(P i) 



j 



\K ( - i 



4ttM 2 
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For the size of the small piece L$ — lnm, damping 
a = 0.05 and 40% polarization degree one gets j rs 
6.7 • 10 7 A/cm 2 using the values of K^ ntr ^ and M for 
cobalt. 
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FIG. 3. Switching phase diagrams for the planar case. 
Regions are marked by letters corresponding the stable 
points. Bistable regions are marked by A/C and A/B. Right: 
time-evolution flow on the projected sphere for bistable A/C 
situation - basin of attraction for C\ is shaded. Upper inset: 
stable points A, B, Ci(I) on the unit sphere. 
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